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MODIFICATIONS TO LANGLEY 0.3-M TCT ADAPTIVE WALL SOFTWARE 
FOR HEAVY GAS TEST MEDIUM 


SUMMARY 


The scheme for two-dimensional wall adaptation with sulfur hexafluoride (SF 0 ) as test gas in the 
NASA Langley Research Center 0.3-m Transonic Cryogenic Tunnel (0.3-m TCT) is presented. 

A unified version of the wall adaptation software has been developed to function in a dual gas 
operation mode (Nitrogen or SF 0 ). The feature of ideal gas calculations for Nitrogen operation is 
retained. For SF 0 operation, real gas properties have been computed using the departure func- 
tion technique. Installation of the software on the 0.3-m TCT ModComp-A computer and 
preliminary validation with nitrogen operation were found to be satisfactory. Further validation 
and improvements to the software, if necessary, will be undertaken at a later date when the 0.3- 
m TCT is ready for operation with SF 0 gas. 

INTRODUCTION 

The Langley 0.3-m TCT 1 is a unique wind tunnel facility for testing two-dimensional airfoil 
models at transonic Mach numbers and flight equivalent Reynolds numbers. The high Reynolds 
number capability is achieved by using cryogenic nitrogen gas as the test medium. The 0.3-m 
since its first proof-of-the concept studies in 1976, has undergone several major changes to 
provide additional testing features. In its present form, the 0.3-m TCT has a 13-mch square test 
section with adaptive wall capability for top and bottom walls and a sophisticated side wall 
boundary- layer removal system. These features make possible testing of airfoil models with a 
greatly reduced level of wall interference. 

The application of the cryogenic concept has been highly successful in simulating the much 
needed high Reynolds number capability at transonic speeds. Many cryogenic wind tunnels have 
been built in the last decade and are being used routinely to test under high Reynolds number 
conditions. Presently, an alternative approach of using a heavy gas in place of the cryogenic 
nitrogen to provide high Reynolds numbers is under consideration. The use of heavy gas per- 
mits operation under ambient temperature conditions making model design, and force and pres- 
sure measurement instrumentation simpler. Also, as demonstrated by Smelt 2 , the power 
requirements will be lower. 
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The present heavy gas feasibility study at Langley uses sulfur hexafluoride (SF 0 ) which is about 
five times heavier than nitrogen/air. However, at a given test Mach number, the increase in 
Reynolds number is about two times due to lower speed of sound ( 30% of air) and lower vis- 
cosity (80%). The use of SF fl also poses much less environmental concerns compared to CFCs. 
However, the effect of the real gas properties of SF 6 on the flow around the model and on the 
shock location at transonic speeds is not well known. The different value of the ratio of specific 
heats and other real gas effects may require corrections to be applied to the test data for valid 
comparison with the nitrogen test data. Therefore, to facilitate a comparative study between 
using cryogenic nitrogen and SF 0 , the 0.3-m TCT is being modified to operate using either gas 
mode. The proposed modification permits a direct tunnel- independent aerodynamic evaluation 
of the two different methods of generating high Reynolds numbers. 

The objective of the present work is to study the aspects of using the 0.3-m TCT adaptive wall 
capability with the SF 0 test gas. Necessary modifications to the adaptive wall testing technique 
in support of the aerodynamic evaluation studies using SF 0 gas have been considered. The adap- 
tive wall software currently in use has been developed specifically for using nitrogen gas as the 
testing medium. For use with SF 0 gas as the test medium, it is necessary to incorporate changes 
to account for the thermodynamic properties of SF 0 and other operational features. The new 
version of the adaptive wall control software has been developed to incorporate the following 
main features. 

a) Easy switching between two test media (Nitrogen & SF 0 ) 

b) No major changes at operator end; thus avoiding operator retraining. 

c) Retain the same program structure and add additional features 
to facilitate easier maintenance of the software. 

This report describes the methodology used in modifying the software for the dual gas mode and 
provides details of the modified software developed during this first phase of the study. Limited 
tests were made to check the performance of the software. Final checking and changes, if 
needed, will be made when the tunnel is ready for operation with SF 0 gas. 

Nomenclature 

, 4 q Constants in expression for constant-pressure ideal heat capacity (Eqn. 2) 
a Speed of sound (meter/second) 
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Constants in equation of state (Eqn. 1) 

Constant in equation of state (Eqn. I) 

Constants in equation of state (Eqn. I) 
Constant-Pressure heat capacity 
Constant- Pressure ideal heat capacity 
Constant- Volume heat capacity 
Constant-Volume ideal heat capacity 
Constants in equation of state (Eqn. 1) 

Enthalpy (JouleAmole) 

Constant in equation of state (Eqn. 1) 

Mach number 
Nitrogen 

Pressure (Newton/meter ) 

Universal gas constant (» 8314.34 joule/(^mole Kelvin) 
Entropy (Joule/(A:moIe Kelvin) 

Sulfur hexafluoride 
Temperature (Kelvin) 

Constant in equation of state (Eqn. 1) 

Velocity (meter/second) 

Volume (meter 3 /A:mole) 

Molecular weight 
Compressibility factor 
Ratio of specific heats (C p /C v ) 

Density (kg/m 3 ) 


Subscripts 

t Stagnation conditions 

v Volume 

1 Static conditions after isentropic expansion 

0.3-m TCT Adaptive Wall Test Section 

The 0.3-m TCT adaptive wall test section (fig. 1) uses a flexible wall concept for wall adaptation. 
It has rigid side walls and flexible stainless steel plates for the ceiling and the floor. The test 
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section has a square cross section with a nominal height of 0.33m (13.0 in.) with the flexible walls 
set to a flat shape. The flexible plates, 1.82m (71.7 in.) long, are anchored at the upstream end 
and are free to move longitudinally in a sliding joint at the downstream end. 

The wall shapes are controlled by stepper motor operated jacks positioned at twenty-one locations 
on the ceiling and floor. The wall positions are monitored by linearly variable differential 
transducers (LVDT) located at each jack station. The wall positions are referenced to .the un- 
deflected straight wall shape. In addition to the wall shapes, the wall streamlining algorithm 
requires the fluid velocity distribution along the wall. This is calculated from measured wall 
pressures at the jack locations from static pressure orifices. The reference test Mach number is 
determined by a static pressure measurement near the upstream anchored end of the flexible 
plates. This station is 0.8m (31.25 in.) upstream of the model mounting turntable center. 

Wall Adaptation Method 

Wall adaptation is the process of adjusting the test section wall shapes to nearly free air stream- 
line shapes. Modern methods developed in the last decade are based on the iterative concept * 
proposed independently by Ferri and Baronty 3 and Sears 4 . The iterative concept treats the flow 
over the model in two parts; the real wind tunnel flow inside the test section and an imaginary 
flow extending to infinity outside the test section walls. For a test section with flexible walls for 
the boundaries, the matching of the flow conditions at the boundaries leads to wall shapes close 
to free air streamline shapes. This concept of streamlining was first introduced by Goodyer 5 to 
define the streamlined conditions for a solid flexible wall wind tunnel. The wind tunnel 
measurements give the flow velocity and direction at the boundaries. One of these quantities 
along with the boundary condition at infinity is sufficient to calculate the imaginary flow. The 
other quantity serves as a compatibility check. Iterative calculation of the real and imaginary 
flows leads to identical flow conditions at the interface boundary. 

There are two approaches to applying the iterative concept of adapting the walls. The first 
method, known as the interface discontinuity method, is particularly suitable to two-dimensional 
airfoil testing. The interface discontinuity method treats the imaginary flow outside the test sec- 
tion ceiling and floor separately. The calculated wall deflections for each wall are then modified 
by adding a portion of the deflection for the opposite wall to represent the coupling effects. 

The second approach 8 considers the coupling effects by using the Cauchy integral formula and 
eliminates the need for separate calculation of the imaginary flows for the ceiling and floor, 
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Further, the concept can be extended to three dimensional test section using Green’s theorem. 


For the 0.3-m TCT, the interface discontinuity approach has largely been used. A limited 
study 7 using the second approach has been made. The purpose of this study was to demonstrate 
the concept of making simultaneous calculation of wall adaptation and residual interference, and 
also to identify some of the shortcomings 8 of the wall adaptation method (WAS-1) then in use at 
the 0.3-m TCT. 

0.3-m TCT Adaptive Wall Software & Present Status 

The 0.3-m TCT adaptive wall software has been developed over a number of years by several 
authors. The software function can be broadly divided into two categories. The first category 
relates to the control of the flexwall hardware, measurement of wall pressures and positions, and 
associated safety features. This portion of the software function will not be affected by the 
heavy gas operation. The second category, critical to the adaptive wall function, is the algorithm 
used to calculate the wall shapes. 

As mentioned earlier, the algorithm used in the 0.3-m TCT is based on the interface discon- 
tinuity approach. The related studies were primarily developed at the University of 
Southampton, England and adapted to the 0.3-m TCT at various stages. Within the concept of 
interface discontinuity, different approaches to calculate the imaginary flow are possible. 
Reference 9 gives the details of a predictive strategy known as WAS-1, first developed at the 
University of Southampton. The WAS-1 method was implemented later at the 0.3-m TCT and 
has been used in a number of airfoil tests. Reference 10 which describes in detail the 0.3-m 
TCT Wall adjustment strategy software, incorporates the WAS- 1 method. 

Despite its wide use in a number of airfoil tests, the experience at Langley in using the WAS-1 
was not encouraging 8 . Hence, during late 1988, the WAS-1 method was replaced with a dif- 
ferent method known as 1FLEX. This major revision is not documented and should be con- 
sidered as an integral part of reference 10. 

Reference 1 1 gives the details of the IFLEX method and a program in HP Basic language. The 
IFLEX method creates an imaginary semi-infinite flowfield which has a streamline originating at 
infinity and has the same shape as the flexible wall. The method represents (fig. 3) the wall 
shapes by a distribution of sources and sinks on a line corresponding to the straight wall shapes. 
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The sources/sink strengths are determined by the local wall slopes. The difference between the 
measured wall pressures and the imaginary flowfield field solution is a measure of new wall 
movements required to adjust the walls to streamline shapes. The procedure is again iterative. 
Since the top and bottom wall shapes are generally different, separate calculations are made for 
each wall. 

The 1FLEX approach overcomes the shortcomings of the original WAS-1 method. Since its in- 
stallation, the streamlining performance has been found to be satisfactory for the airfoils tested. 
However, with a new computer operating system installed during 1991, difficulties were encoun- 
tered with convergence of the wall contours. This necessitated a complete check-out of the 
software before modifications for the SF 6 operation could be incorporated. Details of the 
preliminary check-out program are given in Appendix A. The problem was found to be unre- 
lated to IFLEX performance but required changes to the code to work satisfactorily with the new 
operating system. 


Wall Adaptation for Sulfur Hexafluoride Test Medium 

The principle of wall adaptation as used in the IFLEX approach is essentially based on incom- 
pressible flow theory along with the Prandtl-GIauret rule to take into account the effect of com- 
pressibility. This assumption implies that the flow conditions at the wall be subsonic whereas 
the flow over the model can still have regions of supersonic flows and shock waves. This condi- 
tion is satisfied for most of the transonic testing conditions encountered in the tunnel. 

The requirement of subsonic flow conditions implies that the shape of the free-air streamline 
shapes are independent of the gas properties for a given test Mach number. Therefore, to first 
order, the wall shapes in the fully adapted conditions must be the same for both the nitrogen and 
sulfur hexafluoride operation. However, the approach to the final adapted conditions will be 
different because of the different thermodynamic properties. For a given wall shape and free 
stream Mach number, the measured wall pressures will be different for nitrogen and sulfur 
hexafluoride. Correspondingly, the calculated speed of sound and the fluid velocities at the 
walls will also be different. However, the values of the pressure coefficients remain the same in 
both cases. This feature has been utilized to simulate SF 6 with nitrogen operation. This option 
will be helpful during final validation tests. 
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Results and Discussion 


Properties of Sulfur Hexafluoride 

The change in wall adaptation calculation for sulfur hexafluoride gas arises from the fact that the 
ratio of specific heats is about 1.1 compared to 1.4 for the nitrogen gas. The present adaptive 
wall software for nitrogen operation uses a constant value of 1.4 for all test conditions. This ap- 
proximation is adequate since nitrogen behaves like a perfect gas over the operating conditions of 
interest. 

For SF operation, the following three options are possible for representing the real gas effects in 
wall adaptation calculations. 


a. Use a specified constant value of 7 SF6 

b. Use value of 7g F6 based on stagnation conditions. 

c. Perform complete real gas calculations to calculate 
the thermodynamic properties. 


In contrast to nitrogen, the thermodynamic behavior of SF 6 departs from ideal gas even at fairly 
low pressures. Hence the calculation of thermodynamic properties require solution of the equa- 
tion of state and ideal heat capacity equation 12 : 


P 


RT 

(V - b) 


I 


n. + h.T + c. tw/H-kT/T.) 

(v - by 


C p ° mA o + A 1 T + A z j2 + A s T * + A * j 


( 1 ) 


( 2 ) 


Reference 1 3 gives in detail the solution of the above equations and also the procedure to calcu- 
late properties after isentropic expansion using the departure function technique. In the present 
calculations, the method of reference 13 has been used. Appendix B gives a summary of the 
results obtained in reference 13. 
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Figures 4a-e show the change in thermodynamic properties of SF 6 over the range of pressures (1- 
8 atm.) and temperatures (275-350K) of interest to testing in the 0.3-m TCT. The ratio of 
specific heat, 7, increases from about 1.1 at 1 atmospheric pressure to about 1.13 at 6 atmospheres 
under ambient temperature conditions. Hence, it would be difficult to fix a priori a constant 
value for the ratio of specific heats to simplify adaptive wall calculations. The compressibility 
effects increase rapidly with increase in pressure (fig. 4b). With increase in pressure, the acous- 
tic speed (fig. 4c) decreases at a given temperature. Hence, for a given Mach number the fluid 
velocity will be lower at higher pressures than at lower pressures. The decrease in density with 
increase in temperature is shown in fig. 4d. The specific heat at constant pressure (fig. 4e) 
shows non-linear changes for higher pressures and temperatures below about 275 K. Hence it 
appears that to avoid rapid variations in the properties, it is advisable to work at a stagnation 
temperature of about 300K or higher. 

In addition to these features, the isentropic expansion causes a drop in pressure and temperature. 
The determination of the Mach number after isentropic expansion from stagnation conditions 
requires solution for the entropy equation along with the equation of state and constant pressure 
ideal heat capacity (Appendix B). As mentioned earlier, reference 13 gives a computer program 
to calculate the isentropic expansion. Figure 5 shows the Mach number after isentropic expan- 
sion from an initial stagnation condition of 8 atmospheres pressure and 300K temperature. The 
Mach numbers obtained from real gas calculations are higher than the values using a constant 
value of 7 based stagnation conditions. The differences are smaller at lower Mach numbers. 

The use of heavy gas though requires higher pressures to increase the Reynolds numbers. To ac- 
count properly for departure from ideal conditions, it is necessary to retain real gas calculations 
in the wall adaptation process in the initial studies. 

Simplified calculations based on a fixed value of 7 may be possible after some experience is 
gained. While such a simplified calculation may be adequate from the wall adaptation process, 
the residual errors in the Mach number may be of uncertain magnitude if different types of cal- 
culation are performed for tunnel measurements (ModComp-B) and wall adaptation (ModComp- 
A). Hence, in the present work, complete real gas calculations based on the departure function 
technique have been implemented. Further modifications to provide an option for using a fixed 
value of 7 can be incorporated. 
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Program Implementation 


A Fortran program suitable for implementation on the 0.3-m TCT ModComp-A computer was 
developed based on the method described in reference 13. The program was validated on a HP- 
workstation before installing on the Modcomp computer and integrating with the adaptive wall 
software. The original version of the program presented in reference 13 uses double precision 
calculation to solve the transcendental equations. It was felt that the double precision calcula- 
tion may not be required in a real time environment. Further, double precision is not com- 
patible with the remaining calculations, performed using single precision. 

Table la compares the results of single and double precision calculations on the 0.3-m TCT 
ModComp 16-bit computer for a high pressure test case. The difference between the two results 
is not significant, implying that the single precision calculation is adequate. Table lb and lc 
show examples of single precision calculations for stagnation pressures of 0.6 and 6 atmospheres 
respectively. The results are also compared with results obtained on a HP workstation for valida- 
tion purposes. However, using single precision calculation requires the convergence limits to be 
altered suitably. A typical value of convergence used was 10'® for the ratio of the change in - 
parameter value to the initial value between two iterations. 

Program Integration: 

The real gas program discussed above was integrated with the adaptive wall software described in 
reference 10. The program integration required several changes to many subroutines so that the 
integrated software works in the dual gas mode. Appendix D gives a summary of the various 
programs and subroutines in the modified software. A source listing of the modified and new 
programs added is presented in a separate report. The program is now presently installed on the 
following partitions of 0.3-m TCT ModComp-A computer. 

a. Source files ... Partition ESL 

b. Object files ... Partition EUL 

c. Load files ... Partition ELM 

The start-up and execution procedures remain the same except the load module is assigned to the 
new partition by command //ASS LM ELM. No changes have been made in the OAP system. 


In the present version, the software developed has the following features. 
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a. Functions in both the nitrogen and SF 0 modes. 

b. Choice of the mode, either nitrogen or SF 0 , is made from the digital 
control panel. The values of the parameter IANAL are incremented 
by 10 to represent SF 0 mode of operation (Appendix C) 

c. All the previous features of the control software have been retained. 

Hence, procedures at the operator level remain unaffected. 

d. Automatic or manual selection of scaling factors (Appendix C) 

e. Appropriate message, whether in the nitrogen or SF 0 mode, 
is displayed on the terminal at various stages of program 
execution during testing. 

Appendix C gives a summary of the various settings on the digital control panel for the present 
version of the software. To facilitate installation and checking the various functions of 
the program, a simulation routine has been added. The function of this subroutine is to convert 
the pressures measured during operation with nitrogen to equivalent pressures with SF 0 operation. 
This will be helpful in initial checking and can be retained, if necessary, as an additional feature 
of the software. Suitable changes are necessary after final validation to retain or delete 
this feature. 

Vali dation Tests 

The integrated program was successfully installed and wind-off checks were made to ensure the 
execution of the software in the modified version. A limited number of wind-on checks were 
made to check execution of the flexwall data acquisition and ensure satisfactory execution in the 
nitrogen mode. One of the tests involved a complete wall streamlining cycle with a space shuttle 
model Though the two-dimensional streamlining procedure for a three-dimensional model may 
not be exact, the test was mainly meant to prove the functioning of the software in the nitrogen 
mode.. These initial tests have been useful and have saved considerable development work that 
would have been required later with the heavy gas. Final wind-on proof tests with the SF 0 gas 
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will be conducted later when the 0.3-m TCT is ready for operation with SF # gas. In the present 
form, the real gas calculation portion of the software has optional outputs which will be useful in 
debugging the program. However, the optional printing may slow down the process and will be 
suppressed in the final version. 

Appendix E gives a listing of the program used on the HP workstation to calculate the properties 
of SF fl gas and isentropic expansion. The program can be executed interactively and has provi- 
sions for real gas calculations or for using a specified value of the ratio of specific heats. This 
option can be installed on the 0.3-m TCT ModComp computer, if found necessary, after final 
validation tests. 
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TABLE la 


COMPARISON OF SINGLE & DOUBLE PRECISION CALCULATIONS FOR SF6 
PROPERTIES ON MODCOMP-A (STAGNATION PRESSURE = 6.0 ATM.) 




MODCOMP-A 

MODCOMP-A 



(REAL*4) 

(REAL*8) 

STAGNATION CONDITIONS 




aszssssesssrasasscss: 




d ypurnw / m a ? __ 

_ 

607000.000000 

607000.000000 

T TFMPERATURE K 

= 

305.000000 

305.000000 

V, VOLUME M A 3/KMOLE 

= 

3.905272 

3.905273 

nFkKITY KG/M A 3 . 


37.399185 

37.399183 

2 , COMPRESSIBILITY FACTOR 

= 

0.934786 

0.934786 

CPO, CONS PR IDEAL HEAT CAP. J/KMOLE, K 

= 

98714.062500 

98714.029078 

CVO, CONS VOL IDEAL HEAT CAP. J/KHOLE.K 

- 

90399.718750 

90399.689078 

CP-CV J/KMOLE K 

— 

11429.773437 

11429.783681 

CP, CONS PR HEAT CAP. J/KMOLE, K 

K 

104350.812500 

104350.787244 

CV, CONS VOL HEAT CAR. J/KMOLE ,K 

= 

92921.031250 

92921 .003563 

GAMMA (RATIO CP/CV) 

S 

1.123005 

1.123005 

SOUND SPEED, METERS/SEC 

= 

130.293091 

130.293218 


1SENTR0PIC EXPANSION: 

TstsasssB33sssrss=is= 


P, NEWTON/M A 2 


420000.000000 

420000.000000 

T TEMPERATURE K 

- 

294.678345 

294.678423 

V, VOLUME M A 3/KMOLE 

s 

5.540188 

5.540190 

DENSITY, KG/M A 3 

s 

26.362640 

26.362633 

Z, COMPRESSIBILITY FACTOR .... 

* 

0.949725 

0.949725 

CPO, CONS PR IDEAL HEAT CAP. 

J/KMOLE, K = 

96314.781250 

96314.822471 

CVO, CONS VOL IDEAL HEAT CAP. 

J/KMOLE , K = 

88000.437500 

88000. 482471 

CP*CV, J/KMOLE ,K 


10666.144531 

10666.147160 

CP, CONS PR HEAT CAP. J/KM01E,K = 

100819.187500 

100819.204486 

CV, CONS VOL HEAT CAP. J/KMOLE 

,< = 

90153.031250 

90153.057325 

GAMMA (RATIO CP/CV) 

= 

1.118311 

1.118312 

SOUND SPEED, METERS/SFC 

= 

129.938538 

129.938625 

VEIOCITY (M/SEC) 

r = 

108.851898 

108.851251 

MATH WilMnr P . _ _ 

w 

0.837718 

0.837713 

TEST-SECTION MASS FLOW (KGS/SEC) = 

312.788940 

312.787014 

DYNAMIC PRESSURE (N/M‘2) 

... = 

156181.937500 

156180.082255 

REYNOLDS NO. /METER 

= 

0. 183099E+09 

0.183098E+09 

RATIO PSTAT IC/PTOTAl 

= 

0.691927 

0.691928 

RATIO PSTAT IC/PTOTAL 


0.966158 

0.966159 

RATIO RHO- STAT IC/RHO- TOTAl ... 

= 

0.704899 

0.704899 

RATIO DYN-PR/PTOTAL 

..... = 

0.257301 

0.257298 
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TABU 1 1> 


COMPARISON OF SINGLE PRECISION CALCULATIONS FOR SF6 
(STAGNATION PRESSURE =0.6 ATM.) 


STAGNATION CONDITIONS 

«SSS*3S35S===*a:SB5SS 

P, NEWTON/M '2 

T, TEMPERATURE, K 

V, VOLUME M'3/KMOLE 

DENSITY, KG/M"3 

Z, COMPRESSIBILITY FACTOR 
CPO, CONS PR IDEAL HEAT CAP. J/KM01E,K = 987K. 062500 
CVO, CONS VOL IDEAL HEAT CAP. J/KMOLE, K = 90399.718750 


CP-CV, J/KMOLE, K * 8569.808594 

CP, CONS PR HEAT CAP. J/KMOLE, K « 99208.937500 

CV, CONS VOL HEAT CAP. J/KMOLE, K = 90639.125000 

GAMMA (RATIO CP/CV) * 1.094549 

SOUND SPEED, METERS/SEC = 136.995361 


1SENTROPIC EXPANSION: 

;bci s3rcsx83s:s:z=SB~ 


P, NEUT0N/M*2 E 35000.000000 

T, TEMPERATURE, K = 290.833618 

V, VOLUME M A 3/KMOLE = 68.796326 

DENSITY, KG/M A 3 = 2.122991 

Z, COMPRESSIBILITY FACTOR * 0.995774 

CPO, CONS PR IDEAL HEAT CAP. J/KMOLE ,K = 95394.500000 

CVO, CONS VOL IDEAL HEAT CAP. J/KMOLE, K - 87080.156250 

CP-CV, J/KMOLE ,K = 8491.566406 

CP, CONS PR HEAT CAP. J/KMOLE ,K = 95758.843750 

CV, CONS VOL HEAT CAP. J/KMOLE ,K = 87267.281250 

GAMMA (RATIO CP/CV) = 1.097305 

SOUND SPEED, METERS/SEC s 134.214905 


VFLOCITY (M/SEC) 

MACH NUMBER ■■ 

TEST-SECTION MASS FLOW (KGS/SEC) 

DYNAMIC PRESSURE <N/M A 2) 

REYNOLDS NO. /METER 

RATIO PSTATIC/PTOTAL 

RATIO PSTATIC/PTOTAL 

RATIO RHO-STATIC/RHO TOTAL ...^ 
RATIO DYN-PR/PTOTAI 


= 136.300598 
= 1.015540 
= 31.540771 
= 19720.304687 
= 0. 187199E+08 
= 0.576606 
= 0.953553 
= 0.603478 
3 0.324881 


MODCOMP-A 
(16 bit) 


60700.000000 

305.000000 

41.517105 

3.517923 

0.993775 


PROPERTIES 


HP-345 
(32 Bit) 


60700.000000 

305.000000 

41.517113 

3.517923 

0.993775 

98714.039000 

90399.703000 

8569.810500 

99208.937000 

90639.125000 

1.094549 

136.995460 


35000.000000 

290.833670 

68.796341 

2.122991 

0.995774 

95394.546000 

87080.203000 

8491.565400 

95758.906000 

87267.343000 

1.097305 

134.214980 

136.300330 

1.015538 

31.540716 

19720.232000 

0.187199E+08 

0.576606 

0.953553 

0.603478 

0.324880 
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TABLE lc 


COMPARISON OF SINGLE PRECISION CALCULATIONS FOR SF6 PROPERTIES 


(STAGNATION PRESSURE = 6.0 ATM.) 




MOOCOMP-A 

HP-345 



(16 Bit) 

(32 Bit) 

STAGNATION CONDITIONS 




|S3B2S8tSaB33SrBaCS=: 




P, NEWTON/M '2 

= 

607000.000000 

607000.000000 

T, TEMPERATURE, K 

8 

305.000000 

305.000000 

V, VOLUME M'3/KMOLE 



3.905272 

3.905273 

DENSITY, KG/M'3 


37.399185 

37.399178 

2, COMPRESSIBILITY FACTOR .... 

5 

0.934786 

0.93478 7 

CPO, CONS PR IDEAL HEAT CAP. 

J/KMOLE, K = 

98714.062500 

98714.039000 

CVO, CONS VOL IDEAL HEAT CAP. 

J/KMOLE, K = 

90399.718750 

90399.703000 

CP-CV, J/KMOLE, K 

„ , , T , r , . , « 

11429.773437 

11429.782000 

CP, CONS PR HEAT CAP. J/KMOLE, K * 

104350.812500 

104350.790000 

CV, CONS VOL HEAT CAP. J/KMOLE 

F^ * 

92921.031250 

92921.015000 

GAMMA (RATIO CP/CV) 

T * 

1.123005 

1.123005 

SOUND SPEED, METERS/SEC 

. . . T r . * 

130.293091 

130.293220 

ISENTROPIC EXPANSION: 




l8(E£3=Ef28S3:3£S8=2S 




P, NEWTON/M"2 

s 

420000.000000 

420000.000000 

T, TEMPERATURE, K 

f . . T 8 

294.678345 

294.678300 

V, VOLUME M'3/KMOLE 


5.540188 

5.540190 

DENSITY, KG/M*3 . 


26.362640 

26.362633 

2, COMPRESSIBILITY FACTOR 

= 

0.949725 

0.949725 

CPO, CONS PR IDEAL HEAT CAP. 

J/KMOLE, K = 

96314.781250 

96314.820000 

CVO, CONS VOL IDEAL HEAT CAP. 

J/KMOLE, K = 

88000.437500 

88000.484000 

CP-CV, J/KMOLE, K 

T T . T T - T T * 

10666.144531 

10666.148000 

CP, CONS PR HEAT CAP. J/KMOLE 

,K = 

100819.187500 

100819.210000 

CV, CONS VOL HEAT CAP. J/KMOLE 

f* = 

90153.031250 

90153.062000 

GAMMA (RATIO CP/CV) 


1.118311 

1.118312 

SOUND SPEED, METERS/SEC 

= 

129.938538 

129.938610 

VELOCITY (M/SEC) 

s 

108.851898 

108.851430 

MACH NUMBER 


0.837718 

0.837714 

TEST SECTION MASS FLOW (KGS/SEC) = 

312.788940 

312.787500 

DYNAMIC PRESSURE (N/M'P) 


156181.937500 

156180.590000 

REYNOLDS NO. /METER 

* 

0. 183099E+09 

0.183098E+09 

RATIO PSTAT1 C/PTOTAL 

= 

0.691927 

0.691927 

RATIO PSTATI C/PTOTAL 


0.966158 

0.966159 

RATIO RHO-STATIC/RHO- TOTAL ... 

= 

0.704899 

0.704899 

RATIO DYN-PR/PTOTAl 

” 

0.257301 

0.257299 
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SCHEMATIC OF 0.3-m TCT ADAPTIVE 
WALL TEST SECTION 
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0.3-m Transonic Cryogenic Tunnel Adaptive Wall Test Section 







IMAGINARY FLOW 
(Computed) 




IMAGINARY FLOW 
(Computed) 


Fig. 2 : Real and Imaginary Flow Fields in Wall Adaptation. 
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< N\W 



MODEL 


Fig. 3: Distribution of Sources and Sinks in IFLEX Method. 


19 


Print H GAMMA. pit H RATIO OF SPECIFIC HEATS (SF6) 



Variation of Ratio of Specific Heats for SF 0 with Pressure and Temperature. 




SH 







325 



Print !' CP.pIt i! SPECIF 
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Fig. 4e: Variation of Specific Heat (C p ) for SF 6 with Pressure and Temperature, 




Print H PS-M.pIt j| PS-M (SF6,PT= 8 ATM, TT S 300K) 
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Comparison of Mach number after Isentropic Expansion Using Real Gas 
and Stagnation value of Ratio of Specific Heats. 




APPENDIX A 

Preliminary Check-Out of the 0.3-m TCT Adaptive Wall System 


The purpose of this check-out program was to ensure that the 0.3-m TCT adaptive wall system 
is functioning satisfactorily before the changes for sulfur hexafluoride operation can be imple- 
mented. This effort was necessary in view of the several changes to both the tunnel hardware 
and software systems in the last two years. Initial operation by the 0.3-m TCT staff on a check- 
out model proved unsuccessful in obtaining a successful operation of the adaptive walls and con- 
vergence of wall shapes. Following this effort, the check-out program was initiated using the 
13- inch chord NACA 0012 airfoil model. 

The check-out program consisted of the following: 

a) Primary measurements 

Tunnel total pressure 
Tunnel static pressure 
Tunnel reference pressures 

b) Adaptive wall hardware: 

Wall alignment, LVDT calibrations 

Wall pressure measurements and instrumentation 

c) Control software 

The check-outs for items a-b, and rectifying the observed problems, greatly helped the adaptive 
walls to perform in a predictable manner. Most of the problems were not of a type as to cause 
adverse effects on the wall adaptation process. However, they needed rectification in order to 
achieve the desired test accuracy and trouble free wall movements in a production testing en- 
vironment. One of the major source of errors which affected blockage corrections was a leakage 
in the total pressure measurement tubing. The wall pressure measurements were also checked 
using both tunnel static and vacuum as reference pressures. 

After these changes, it was noted that the walls moved correctly in the first iteration, but never 
seemed to converge in subsequent iterations. The solutions for wall shapes showed continuous 
driving of the walls towards free jet conditions. It was difficult to attribute this behavior to any 
known past trends or experience. No obvious errors or doubts about measurement systems were 
present at this stage. 
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More detailed off-line checks of the wall adjustment software on a HP workstation showed that 
the wall shape solutions were satisfactory and the interference levels low enough to be acceptable. 
This enigma was resolved by tracing the cause of the error to a limitation of the new operating 
system recently installed on the ModComp computer in the 0.3-m TCT. The new system inter- 
preted multiplication of an integer constant and a real variable differently in the subroutine 
IFLEX which calculates new wall shapes for each iteration. This limitation resulted in calcula- 
tion of wall movements in only one direction by ignoring the sign of the real variable. The same 
subroutine had worked satisfactorily under the old operating system. This spurious calculation 
resulted in unidirectional movement of the walls, increasing the level of residual interference 
with additional iterations showing no sign of convergence of wall shapes. 

This limitation was overcome by defining the integer constant as a real number in the IFLEX 
subroutine. With this change, the adaptive wall software worked satisfactorily. The successful 
operation was demonstrated on a 13- inch chord NACA-0012 airfoil model at two angles of at- 
tack. Following this, another airfoil model has been tested successfully by the 0.3-m TCT staff 
over a range of Mach numbers and angles of attack. The observed limitation of the new operat- 
ing system requires caution in developing new programs until the system is modified. 

Improvements to software 

During the course of the check-out program, a new procedure for setting the scaling factors for 
top and bottom walls was suggested. The new procedure provides a choice for choosing scaling 
factors either automatically or manually through the digital constants panel. This change 
eliminates the need for maintaining different programs for these tasks. Often, while testing un- 
der hard to adapt conditions, it is desirable to use the manual feature to set lower values of the 
scaling parameters to improve convergence. Mr. Eric Walberg implemented the new procedure 
in the adaptive wall software. 
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APPENDIX B 

Isentropic Expansion of Sulfur Hexafluoride (SF6) 

• f 12 is 

The equation of state and the ideal constant-pressure heat capacity for SFg are given by 


P 


i=6 

RT 

(V - b) 

i=2 



a. + /».T + c { exit- kT/T e ) 
(V - h y 


<7 


= a +A.T + A.;r + 

o l * 




A 4 T 


The various constants equations A I and A2 are given by 


B1 


B2 


h 

sac 

0.047812001 



k 

= 

6.8830220 



T 

C 


318.8 



a 2 

a- 


- 1.064506759 
1.284952625 

X 

X 

It) 6 

I0 5 

3 

a A 

zz 

- 7.338851897 

X 

10 s 

»; 

s 

- 3.25590376 




m 

1.170089157 

X 

10 s 


= 

-1. 040549 130 

X 

10 2 

h 4 


0.0 



*5 

= 

7.271 143381 

X 

10 



- 5.067990044 

X 

I0 7 

l 


8.783983640 

X 

I0 6 

i\ 

— 

0.0 



r r» 


- 2.0485.35876 

X 

I0 4 

A a 

= 

- 15748.49323 



O 

A i 

A 2 


575.7699892 



= 

_ 0.749044588 



< 

3S 

3.538653215 

X 

10' 

= 

- 1.402388775 

X 

It) 8 


I'or given values of pressure and temperature, the equation ol state is solved to calculate the 
volume V. The constant pressure ideal heat capacity depends only on the temperature. The 
corresponding constant-volume ideal heat capacity C y ° is given by 


C 


O 


V 


C " - R 

i> 


B3 
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where R is the universal gas constant. 


The corresponding expressions for the heat capacities, C p and C y for the non-ideal behavior of 
the gas are given by 

Cp - C y - -T (3P/3T) 2 /(3P/aV) B4 

C y - C v ° + T {d z P/dT\dy B5 

For the real gas expansion from stagnation conditions at temperature T t and volume V t , to condi- 
tions at temperature and volume V v the change in enthalpy (H) and entropy (S) are given by 


pTt fVt fV t 

H(T t ,V t ) - H(T r Vp - j T C p ° JT - j JP t - RT t /V t ) d\ t + T t O/3T t )J (P t -RT t /V t ) d\ 

+ f V1 (P x - RTj/Vj) d\ l - T 1 (3/aT 1 )f V1 (P -RT./V.) dV. 

j OO J oo 


+ (P t V t -RT t ) + (PjVj-RTj) 


B6 


S( 


T t ,V t ) - S(Tj,Vj) rf Tt (Cp°/T) dT - R log( T t /V t ) + (3/3T t ) f V ‘(P t -RT t /V t ) dV t 
T1 J oo 


- R /oi'fTj/Vj) + (d/dT i f V1 (P 1 -RT 1 /V 1 ) rfVj 

J oo 

B7 

The equations A6 and A7 are derived using the departure function technique defined by the dif- 
ference in gas properties between the real and ideal states. For the particular case of isentropic 
expansion, the change in entropy between the initial and final conditions is zero. The tempera- 
ture after isentropic expansion is determined by solving equations for state, heat capacity and the 
entropy. Because of the transcendental nature, the equations are solved iteratively. Knowing 
the change in enthalpy, the fluid velocity U after isentropic expansion is calculated from. 


U = ( 2 (H(T t ,V t ) - HCT^jH/WI 1 '* 


B8 


where VV is the molecular weight. 
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APPENDIX C 


DIGITAL CONTROL PANEL SETTINGS 


Parameter 

Symbol 

N 2 

SF 6 

Remarks 

Analysis 

1ANAL 

0 

10 

Airfoil Tests 



1 

11 

Empty TS Streamlining 



2 

12 

Re-analysis (Airfoil Tests) 



3 

13 

Development Streamlining 

Scaling 

TWSF & 

0.0 

0.0 

Automatic Scaling 

Factor 

BWSF 



Iteration el 0.80 

Iteration =2 0.70 

Iteration>=3 0.60 



>0.0 

>0.0 

Set Values 

For all Iterations>=l 

Coupling 

TWC'PLF 

0.35 

0.35 

Recommended Value 

Factor 

BWCPI.F 

— 


(Do not change) 
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APPENDIX D 


SUMMARY OF ADAPTIVE WALL SOFTWARE MODIFICATIONS 


Name 

Status 

Changes 

Remarks 

FLXWAS* 

E 

Yes 

Statements & calls related 
to IANAL parameter 

INITWA + 

E 

No 


STWAl.L 

E 

No 


WALCAL + 

E 

No 


STAR + 

E 

No 


WALDAT 

E 

Yes 

Statements & calls related 
to IANAL parameter 

WAS 

E 

Yes 

Statements & calls related 
to IANAL parameter 

SUME 

F. 

Yes 

Statements & calls related 
to IANAL parameter 

OUT 

E 

Yes 

Statements & calls related 
to IANAL parameter 

IFLEX 

F.,M 

No 

Not documented. Further 


Modifications Nov ’91. 


ERROR + 

SMOOTH 

FILESORT 

GETDATA 

LIFT 

BLOCK 

E 

E 

E 

E 

E 

E 

No 

No 

No 

No 

No 

No 

Sorts array in ascending order 
Gets wall data from Ref. Table 

OAPCM 

E 

No 

Common Definition (System) 

FLXCOM 

E 

No 

Common Definition (Flexwall) 

WASCOM 

E 

No 

Common Definition (WAS) 

FLXTYP 

E 

No 

Type Declaration (Flexwall) 

FWPTYP 

E 

No 

Type Declaration (WASPAR) 

WASTYP 

E 

No 

Type Declaration (WAS) 

FWPEQU 

E 

No 

Equivalence (WASPAR) 

FMASSIGN 

E 

No 

Reference Table 


continued. 
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APPENDIX D 


SUMMARY OF ADAPTIVE WALL SOFTWARE MODIFICATIONS 


Name 

Status 

Changes Remarks 

STARSF 

New 

Boundary-layer calculation 
for SF 0 gas 

SIMUL 

New 

Simulation of SF 0 with N 2 . 
(For initial checking) 

SF6 

New 

Mach Number, Velocity and 
Density Calculation (SF„) 
Isentropic Expansion (SF 6 ) 

SISEN 

New 

SVOL 

New 

Solves equation of State (SF g ) 

SSPHEAT 

New 

C , Specific heat at Constant 
Pressure Calculation. (SF g ) 

SINTCP 

New 

Calculates Integral of C from 
Temperature T x to T t (ST 6 ) 

SENRCP 

New 

Integral for Entropy change 
calculation (SF 0 ) 

SRCONS 

New 

Calculation of constants (SF 6 ). 


* Main Program; E: Existing; M: Modified 
+ I/O changes only to compile on HP-Workstation. 


Concluded. 


Note: In addition to the programs listed in the table, the adaptive wall software 

contains/uses several other system calls and procedures for dedicated tasks. 
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APPENDIX E 

PROGRAM FOR CALCULATING SF fl PROPERTIES 


PROGRAM MAIN 

IMPLICIT REAL* 8 (A-H,0-Z) 
DIMENSION PSTAT(22) , EM(22) 


DENS (22 ) , VELFPS (22) , IERR(22) 


PRINT *, 'ENTER TOTAL PRESSURE IN PSIA; ' 
READ * , PTOTAL 

PRINT *, 'ENTER TEMP IN K ; 

READ * , TTOTAL 

PRINT *, 'NUMBER OF STATIC PRESSURES;' 
READ * , NPORTS 


PRINT *, 'ENTER STATIC PRESSURE IN PSIA; ' 

PRINT *, 'INPUT NEGATIVE VALUE TO CALCULATE FOR STAG ONLY) ; 

DO 100 1=1, NPORTS 

PRINT *, I 

READ * , PSTAT ( I ) 

100 CONTINUE 

PRINT *, 'ENTER PRINT PARAMETER, 0,1 OR 2; ' 

READ *, IPRN 

PRINT *, 'ENTER TYPE OF CALCULATION;' 

PRINT *, '0 FOR REAL GAS CALCULATIONS' 

PRINT *, '1 FOR IDEAL GAS (GAMMA STAG VALUE)' 

PRINT *, '2 FOR IDEAL GAS (SPECIFIED GAMMA VALUE)' 

READ *, ICODE 

IF ( ICODE. EQ . 2 ) READ * , SPGAMA 


CALL SF6 ( PTOTAL , TTOTAL , NPORTS , PSTAT , ICODE , I PRN , 
+SPGAMA, EM, DENS, VELFPS, I ERR) 

PRINT *, EM, DENS, VELFPS 
DO 200 1=1, NPORTS 

WRITE (LO, 150) I , PSTAT ( I ) , EM(I) ,DENS(I) ,VELFPS(I) 
150 FORMAT ( 5x, 12 , 4F15 . 6 ) 

200 CONTINUE 
STOP 
END 


C 


C 

c 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE SF6 (PTOTAL, TTOTAL , NPORTS , PSTAT, ICODE, IPRN, 
+SPGAMA, EM, DENS, VELFPS, I ERR) 

IMPLICIT REAL *8 (A-H.O-Z) 

DIMENSION PSTAT (22) , EM(22) , DENS(22) ,VELFPS(22) , IERR(22) 
PURPOSE: 

To'cALCULATE LOCAL MACH NUMBERS, DENSITY AND VELOCITIES 
AT WALL MEASUREMENT POINTS ASSUMING ISENTROPIC EXPANSION. 

METHOD: 

TWO OPTIONS ARE AVAILABLE. A) REAL GAS CALCULATIONS FOR 
EACH PORT. B) ISENTROPIC EXPANSION BASED ON GAMMA VALUE 
FOR STAGNATION CONDITIONS PERFORMED. 

INPUT PARAMETERS : 

PTOTAL : TOTAL PRESSURE (PSIA) 

TTOTAL : TOTAL TEMP (K) 

NPORTS : NUMBER OF PORTS FOR CALCULATION 

PSTAT (NPORTS) : STATIC PRESSURE (PSIA) NPORTS TIMES 

. INPUT -VE VALUE FOR STAGNATION CALC ONLY 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ICODE 


IPRN 


SPGAMA 


CALCULATION CODE 

=0 REAL GAS CALCULATIONS 

=1 EXPANSION BASED ON STAGNATION GAMMA 

=2 EXPANSION BASED ON SPECIFIED GAMMA 

PRINT PARAMETER 

=0, NO PRINT OUTOUT 

=1, OUTPUT FINAL PARAMETERS { ICODE=0 ) 

=2, LONG OUTPUT (FOR CHECK) ( ICODE=0 ) 

RATIO OF SPECIFIC HEATS (ONLY FOR ICODE=2) 


RETURN PARAMETERS: 


EM(NPORTS) 

DENS (NPORTS) 
VELFPS ( NPORTS ) 
I ERR (NPORTS) 


/ISENX/ 

/TOTL/ 


LOCAL MACH NUMBERS 

DENSITY ( SLUGS/FT A 3 ) 

VELOCITY (FT/SEC) 

ERROR CODE 
=0 NORMAL RETURN 

=1 PSTAT TOO LOW FOR ISEN CALCULATIONS 
PROPERTIES AFTER ISEN EXPANSION 
PROPERTIES FOR STAG. CONDITIONS 


SUBROUTINES CALLED: 


c 

c 

VOL: 


c 


RCONS 

c 


SPHEAT 

c 

ISEN: 


c 


VOL: 

c 


RCONS 

c 


SPHEAT 

c 


RCONS 

c 


ENRCP 

c 

c 


JNTCP 


COMMON /ISENX/ VEL , FMACH , FMASS , DYNPR , REY , 

+ PS BY PT, T1BYTT, VTBYV1 , DYNBYPT, 

+ T1 , VI , Zl , GAM1 , AS1 , DENI 

COMMON /TOTL/ ZT , GAMT , AST , DENT 
DATA LO/0LO/ 

LO=6 

DO 400 1=1, NPORTS 
IERR( I) =0 

PT=6894.7572D+00*PTOTAL 
PS=6894 . 7572D+00*PSTAT( I) 

TT=TTOTAL 


IF ( IPRN.NE . 0 ) WRITE(LO, 305) 

305 FORMAT(/10X, 'STAGNATION CONDITIONS' , /10X, 


IF(I.EQ.l) CALL VOL ( PT, TT, VT, ZT, GAMT, AST, DENT, IPRN) 
PSCHEK=PS/PT 

IF ( PSCHEK . LT .0.1) GOTO 330 


.AND. ICODE. EQ.0) WRITE(LO, 310) 

.AND. ICODE. EQ.l) WRITE(LO, 311) GAMT 
.AND. ICODE. EQ. 2) WRITE (LO, 312 ) SPGAMA 
ISENTROPIC EXPANSION (ICODE=0; REAL GAS):' 


IF( IPRN.NE. 0 
IF (IPRN.NE. 0 
IF (IPRN.NE. 0 

310 FORMAT ( /10X, 

+ ,/10X, 

311 FORMAT(/10X, 'ISENTROPIC EXPANSION (ICODE=l; STAG GAMA= ' , F8 . 4 , ' ) ' 

+ , /iox, ' ================================*================='; 
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on on no onn 


312 FORMAT(/10X, ' ISENTROPIC EXPANSION (IC0DE=2; SPEC GAMA= ' , F8 . 4 , ' ) ' ^ 


REAL GAS CALCULATIONS FOR EXPANSION 
IF { ICODE . EQ . 1 .OR. IC0DE.EQ.2) GO TO 315 
IF ( ICODE . EQ . 0 ) CALL ISEN ( PT, TT, VT, PS, GAMT, IPRN) 

GOTO 320 

ISENTRPIC EXP BASED ON STAGNATION CONDITIONS 
315 IF ( ICODE . EQ . 2 ) GOTO 318 
G1=(GAMT-1 .0) /2 .0 
G2=(GAMT-1.0) /GAMT 

FMACH=SQRT ( ( (PT/PS) ** (G2) -1.0 ) /Gl ) 

03=1.0/ (GAMT-1.0) 

G =1.0 + G2 * FMACH * FMACH 
DEN1=DENT/ (G**G3) 

AS1 =AST /SQRT(G) 

VEL =FMACH*AS1 
GOTO 320 

ISENTROPIC EXPANSION USINF SPECIFIED GAMMA VALUE 
318 Gl= (SPGAMA-1 . 0 ) /2 . 0 

G2= (SPGAMA-1 . 0 ) / SPGAMA 

FMACH=SQRT ( ( (PT/PS) ** (G2) -1.0 ) /Gl ) 

G3=l .0/ (SPGAMA-1 .0) 

G =1.0 + G2* FMACH *FMACH 
DEN1=DENT/ (G**G3) 

AS1 =AST /SQRT (G) 

VEL =FMACH*AS1 

STORE VALUES 
320 EM ( I ) =FMACH 

DENS ( I ) =DENl/515 . 379D+00 
VELFPS(I) =VEL/0 .3048 
GOTO 400 
330 I ERR ( I ) = 1 
400 CONTINUE 
RETURN 
END 
C 

SUBROUTINE ISEN ( PT, TT, VT, PS, GAMT, IPRINT) 

IMPLICIT REAL *8 (A-H.O-7.) 

DIMENSION AT ( 5 ) ,DAT(5) , DDAT ( 5 ) ,A1(5) ,DA1(5) ,DDA1(5) 

C 

C PURPOSE : 

C ======== 

C TO CALCULATE ISENTROPIC EXPANSION FOR REAL GAS FROM 

C GIVEN STAGNATION PRESSURE S, TEMPERATURE TO A SPECIFIED 

C EXPANSION PRESSURE. 

C 

C ENTRY PARAMETERS: 

C ================= 

C PT: STAGNATION PRESSURE (N/M A 2) 

C TT: STAGNATION TEMPERATURE (K) 

C PS: STATIC PRESSURE AFTER EXPANSION (N/M A 2) 

C VT: STAGNATION VOLUME 

C GAMT: RATIO OF SPECIFIC HEATS FOT PT & VT 

C IPRINT: OUTPUT PARAMETER 

C =0 NO OUTPUT TO DISPLAY 

C =1 RESULTS OUTPUT TO DISPLAY 


1 

I 


35 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 


c 


RETURN PARAMETERS ( IN / ISENX/ ) 


VEL : 
FMACH : 
FMASS : 
DYNPR : 
REY s 
PSBYPT : 
T1BYTT : 
VTBYVl : 
DYNBYPT: 
T1 : 
VI : 
Zl : 
GAM1 : 
AS1 : 
DENI : 


VELOCITY (M/SEC) 

MACH NUMBER 

TEST-SECTION MASS FLOW (KGS/SEC) 

DYNAMIC PRESSURE (N/M A 2) 

REYNOLDS NO. /METER 
RATIO PSTATIC/PTOTAL 
RATIO TSTATIC/TTOTAL 
RATIO RHO- STATIC /RHO -TOTAL 
RATIO DYN-PR/PTOTAL 

STATIC TEMPERATURE AFTER EXPANSION (K) 

STATIC VOLUME AFTER EXPANSION (N/M A 2) 
COMPRESSIBILITY FACTOR 

RATIO OF SPECIFIC HEAT AFTER EXPANSION 
SPEED OF SOUND FOR STATIC CONDITIONS (M/SEC) 
DENSITY <KG/M A 3) 


SUBROUTINES CALLED 


VOL: 

RCONS 
SPHEAT 
RCONS : 

ENRCP: 

INTCP: 

COMMON /ISENX/ VEL , FMACH , FMASS , DYNPR , REY , 

+ PSBYPT, T1BYTT, VTBYVl f DYNBYPT , 

+ T1 , VI , Zl , GAMl , AS1 , DENI 

LO=6 

CALL RCONS ( TT , AT , DAT , DDAT ) 

XT = 1 . 0D+00/ (VT-0 . 047812001D+00 ) 

ISENTROPIC EXPANSION 
Gl=(GAMT-l . 0D+00 ) /GAMT 
T1=TT* ( ( PS/ PT) * *G1 ) 

G3=2 . 0D+00/ (GAMT- 1 . 0D+00 ) 

FM= (PT/PS) **G1 - l.ODfOO 
FMT=DSQRT(G3*FM) 

TLO=Tl 

TUP=0 ,99999999999D-»*00*TT 
100 VI =8 3 14 . 34D+00*T1/PS 

IF ( I PRINT . EQ . 2 ) PRINT *, T1 , VI 
CALL ENRCP (TT f Tl # SI) 

150 S2 = 8314-34D+00*DLOG((TT*V1)/(VT*T1)) 

XI = 1 . 0D+00/ (Vl-0. 047812001 D+ 00) 

CALL RCONS (T1 , Al , DA 3 # DDAl ) 


4 -DAT ( 1 ) *DLOG (1 . 00D400 + XT*0 . 047812001D+00) 

4 +DA1 ( 1 ) *DLOG ( 1 . 00D+00 + Xl *0 . 0478120 01D+00 ) 

+ - (DAT ( 2 ) *XT - DAI ( 2 ) *X1 ) 

+ -(DAT(3) *XT*XT - DAI ( 3 ) *Xl*Xl ) /2 . 0D+00 

4 - (DAT(4) *XT*XT*XT - DAI (4) *X1*X1*X1) /3 . 0D+00 

4 - (DAT(5) *XT*XT*XT*XT - DAI (5) *X1*X1*X1*X1) / 4 . 0D+05 

IF(IPRINT.EQ.2) PRINT *, 'SI, S2, S3 ' 

IF ( IPRINT . EQ . 2 ) PRINT *, SI, S2, S3 
G = SI - S2 + S3 


36 



c 

c 

c 

c 


c 

c 


c 


DPV = DAI (5) * (Xl**5. OD+OO) + DAI ( 4 ) * (Xl**4 . OD+OO ) 

+ +DA1 (3) *(X1**3 .OD+OO) 

+ +DA1 (2) * (Xl**2 .OD+OO) + DAI ( 1 ) * (XI ) 

DGV = -DPV 

IF(IPRINT.EQ.2) PRINT *, DPV, DGV 
V2 = VI - G/DGV 
IF (I PRINT. EQ. 2) PRINT *, VI, V2 
IF (DABS (V2-V1) .LT. 1 .D-ll) GOTO 200 change date 7/13/92 
IF (DABS (V2-V1) .LT.l. D-ll -OR. V2.LT.VT) GOTO 200 
VI = V2 
GOTO 150 

200 PS1= A1 (5) * (Xl**5. OD+OO) + Al (4) * (Xl**4 .0D+00) 

+ +A1(3)*(X1**3. OD+OO) 

+ +A1(2)*(X1**2. OD+OO) + Al ( 1) * (XI ) 

IF(IPRINT.EQ.2) PRINT *, PS1.PS 

ACCURACY CHANGED 7/17/1992 TO 10 PASCALS TO IMPROVE 
CONVERGENCE AT LOW PRESSURES 
IF (DABS ( PS1-PS) .LT. 0 . 05D+00) GOTO 300 
IF(DABS(PSl-PS) .LT. 10 .OD+OO) GOTO 300 

IF ( ( PS-PS1 ) . LT .0.0) GO TO 250 
TL0=T1 

T1 = (TLO+TUP) /2 .OD+OO 
GOTO 100 
250 TUP=T1 

T1 = (TLO+TUP) /2 .OD+OO 
GOTO 100 
300 CONTINUE 

CALL VOL ( PS, T1 , VI , Z1 , GAM1 , AS1 , DENI , I PRINT) 

ENTHALPY CALCULATION 
CALL INTCP (TT, T1 , SI ) 

52 = -AT ( 1 ) *DLOG ( 1 . 0 0D+00 + XT*0 . 047812001D+00 ) 

+ -AT (2) *XT 

+ -AT(3) *XT*XT/2 .0D+00 

+ -AT(4) *XT*XT*XT/3 . OD+OO 

+ -AT(5) *XT*XT*XT*XT/4 . 0D+00 

PRINT * SI S2 

53 = -Al ( 1 ) *DLOG ( 1 . 0 OD+OO + X1*0 . 047812001D+00) 

+ -Al (2) *X1 

+ -Al(3)*Xl*Xl/2. OD+OO 

+ -Al (4) *Xl*Xl*Xl/3 .OD+OO 

+ -Al(5)*Xl*Xl*Xl*Xl/4. OD+OO 

PRINT *, S3 

54 = -DAT ( 1 ) *DLOG ( 1 . 00D+00 + XT*0 . 047812001D+00 ) 

+ -DAT ( 2 ) *XT 

+ -DAT ( 3 ) *XT*XT/2 . OD+OO 

+ -DAT (4) *XT*XT*XT/3 .OD+OO 

+ -DAT ( 5 ) *XT*XT*XT*XT/ 4 . OD+OO 

54 = TT*S4 
PRINT *, S4 

55 = -DAI (1) *DLOG(l .OOD+OO + X1*0 . 047812001D+00) 

+ -DAI ( 2 ) *X1 

+ -DAI (3)*Xl*Xl/2. OD+OO 

+ -DAI (4) *Xl*Xl*Xl/3 .OD+OO 

+ -DAI (5) *Xl*Xl*Xl*Xl/4 . OD+OO 


S5 = TV 
PRINT ’ 


S5 
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o o on 


S6 = PT*VT - PS*V1 + 8314 . 34D+00* (T1 -TT) 

C PRINT *, S6 

DH = SI - S2 + S3 + S4 - S5 + S6 
IF ( IPRINT . EQ . 2 ) PRINT *, DH 
VEL= DSQRT(2 . 0D+00 *DH/ 146 . 054D+00) 

Gl= (GAMl-i . OD+OO) /GAM1 
G3=2 . OD+OO/ (GAM1-1 . 0D+00 ) 

FM= ( PT/PS ) * *Gl - 1. OD+OO 
FM1=DSQRT(G3*FM) 

FMACH= VEL/AS1 

FMASS=VEL*DEN1*0 . 1090D+00 

DYNPR=0 . 5D+ 0 0 *DEN1 *VEL*VEL 

REY =DEN1*VEL/ (Tl*5 . 59D-08 - 8.0D-07) 

PSBYPT=PS/PT 
T1BYTT=T1/TT 
VTBYVl =VT/ V 1 
DYNBYPT=DYNPR/ PT 

IF( IPRINT. EQ.l .OR. IPRINT. EQ. 2) WRITE (LO, 3 10 ) 

+VEL, FMACH, FMASS, DENI , DYNPR , REY, PSBYPT, T1BYTT, VTBYVl, DYNBYPT 


310 FORMAT (/10X, 'VELOCITY (M/SEC) =',F20.6, 

+ /I OX, 'MACH NUMBER =',F20.6, 

+ /10X, 'TEST-SECTION MASS FLOW (KGS/SEC) =',F20.6, 

+ /10X, 'DENSITY, KG/M A 3 =',F20.6, 

+ /I OX, 'DYNAMIC PRESSURE (N/M~2) =',F20.6, 

+ /I OX, 'REYNOLDS NO. /METER =',E20.6, 

+ /10X, 'RATIO PSTATIC/PTOTAL =',F20.6, 

+ /10X, 'RATIO TSTATIC/TTOTAL =',F20.6, 

*■ /10X, 'RATIO RHO-STATIC/RHO-TOTAL =',F20.6, 

+ /10X, 'RATIO DYN-PR/PTOTAL =',F20.6) 

C 


IF ( IPRINT. EQ. 2) PRINT *, ' SI , S2 , S3 , S4 , S5, S6 , DH/146 . 054D+00 ' 

IF ( IPRINT. EQ. 2) PRINT *, Si , S2, S3 , S4 , S5, S6, DH/146 . 054D+00 
IF ( IPRINT. EQ. 2) PRINT *, 'VEL, FMACH, FMT, FM1 ' 

IF ( IPRINT. EQ. 2) PRINT *, VEL, FMACH, FMT, FM1 

RETURN 

END 

SUBROUTINE VOL(P,T,V, Z, GAM, AS, DEN, IPRINT) 

C P: IN PASCALS (N/M* *2 ) ; T: IN KELVIN; V: M**3/KMOLE (OUTPUT) 

IMPLICIT REAL *8 (A-H,0-Z) 

DIMENSION A ( 5 ) , DA ( 5 ) , DDA ( 5 ) 

LO=6 

CALL RCONS (T, A, DA, DDA) 

X = P/ (8314.34D+00*T) 

100 F = A ( 5 ) * (X**5. 0D+00) + A(4) * (X**4 . OD+OO) + A{3) *(X**3 .0D+00) 
+ +A(2)*(X**2. 0D+00) + A (1 ) * (X) - P 

FD= 5.0D+00*A(5)*(X**4.0D+00) + 4 . 0D+00*A(4) * (X**3 . OD+OO ) 

+ +3 . 0D+00*A (3)*(X**2. OD+OO) + 2 . OD+00*A(2) *X + A(l) 

Xl= X - F/FD 

IF (DABS(Xl-X) .LE. 1.0D-10) GOTO 200 
X = XI 
GOTO 100 
200 CONTINUE 

CALCULATE VOLUME t* COMPRESSIBILITY FACTOR 
V = 0 .047812001D+00 + (1.0D+00/X) 

Z = (P*V) / (8314 . 34D+00 *T) 

CALCULATE CP-CV 

DPV = -8 . 0D+00*A(5) * (X**6 . 0D+00) - 4 . 0D+00*A(4) * (X**5 . 0D+00 ) 
-3 . 0D+00*A ( 3 ) * (X**4. OD+OO) 

-2 .0D+00*A(2) * (X* * 3 .OD+OO) - A{1) * (X**2 .OD+OO) 
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on on 


DPT = DA (5) * (X**5.0D+00) + DA(4) * (X**4 .00+00) 
+ +DA (3)*(X**3.0D+00) 

+ +DA(2) * (X**2 . 0D+00) + DA ( 1 ) * ( X ) 

CPMCV= -T* (DPT*DPT) /DPV 

DINT IS INTEGRAL OF DDPDTT W.R.T V 
DINT = - (DDA(5) /4 . 0D+00) * (X**4 .OD+OO) 

+ - (DDA(3) /2 . OD+OO) * (X**2 . OD+OO) 

+ - (DDA (2)/1.0D+00)*(X) 

CALL SPHEAT (T, CPO ) 

CVO = CPO - 8314.34D+00 
CV = CVO + T*DINT 
IF ( IPRINT . EQ . 2 ) PRINT *, CPO, CVO 
GAM = 1. OD+OO + (CPMCV/CV) 

CP » CPMCV + CV 


SPEED OF SOUND 

AS = DSQRT{-GAM*V*V*DPV/146 . 054D+00) 

DEN = 146 . 054D+00/V 

IF ( IPRINT. EQ . 1 .OR. IPRINT. EQ. 2) WRITE(LO, 300) 

+P, T, V, DEN, Z, CPO , CVO , CPMCV, CP, CV, GAM, AS 

300 FORMAT ( 10X,'P, NEWTON/M A 2 

4 /10X,'T, TEMPERATURE, K 

4 /10X, 'V, VOLUME M~3/KMOLE 

4 /I OX, 'DENSITY, KG/M A 3 

4 /10X, 'Z, COMPRESSIBILITY FACTOR 

+ /10X, 'CPO, CONS PR IDEAL HEAT CAP. J/KMOLE.K 

+ /10X, 'CVO, CONS VOL IDEAL HEAT CAP. J/KMOLE.K 

4 /10X, 'CP-CV, J/KMOLE.K 

4 / 10X, 'CP, CONS PR HEAT CAP. J/KMOLE.K 

4 / 10X, 'CV, CONS VOL HEAT CAP. J/KMOLE.K 

4 /10X, 'GAMMA (RATIO CP/CV) 

+ /10X, 'SOUND SPEED, METERS/SEC 

RETURN 

END 1 


=' ,F20.6, 
= ' , F20 . 6, 
= ' , F2 0 . 6 , 
=' ,F20.6, 
=' , F20.6, 
=' ,F20.6, 
= ' , F20 . 6 , 
=' , F20.6, 
=' ,F20.6, 
= ' , F20 . 6 , 
=' , F20.6, 
= ' , F20 . 6) 


C 


C 

C 

C 

C 

C 


SUBROUTINE SPHEAT(T,CP) 

IMPLICIT REAL *8 (A-H.O-Z) 

PURPOSE: TO CALCULATE CONSTANT PRESSURE IDEAL HEAT CAPACITY 

ENTRY : T, TEMPERATURE (K) 

RETURN : CP, SPECIFIC HEAT J/(KMOLE.K) 


CP = -15748. 49323D400 + 

4 

4 + 

4 

RETURN 

END 


575.7699892D400 
0 .749044588D+00 
3 . 53 86532 15D-04 
1 .402388775D+08 


* T 

*(T**2. OD+OO) 
*(T**3. OD+OO) 
*(T**(-2. OD+OO) ) 


SUBROUTINE INTCP (TT, T1 , SI ) 

IMPLICIT REAL *8 (A-H.O-Z) 

SI = -15748 .49323D+00* (TT-T1) 

+575.7699892D+00* (TT*TT-T1*T1 ) /2. OD+OO 
-0.749044588D+00* (TT*TT*TT-T1*T1*T1) /3 .OD+OO 
+3 . 538653215D-04* (TT*TT*TT*TT-T1*T1*T1*T1) /4 .OD+OO 
-1. 402388775D+08*(-l. 0D+00/TT +1.0D+00/T1) 

RETURN 

END 
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c 


c 


c 


c 


c 


SUBROUTINE ENRCP (TT, T1 , SI ) 

IMPLICIT REAL *8 (A-H,0-Z) 

SI = - 15748 . 49323D+00*DLOG (TT/T1 ) 

+ 575.7699892D+00*(TT -Tl) 

- 0.749044588D+00* (TT*TT -Tl*Tl) /2 . 0D+00 

+ 3 . 538653215D-04* (TT*TT*TT -Tl*Tl*Tl) /3 .0D+00 
+ 1. 4023887 75D+08*(0. 5D+00/TT**2 . OD+OO) 

- 1. 4023887 75D+08*(0. 5D+00/Tl**2 .0D+00) 

PRINT *, TT, Tl, SI 

RETURN 

END 


SUBROUTINE RCONS (T, A, DA, DDA) 

IMPLICIT REAL *8 (A-H.O-Z) 

DIMENSION A ( 5 ) , DA ( 5 ) , DDA ( 5 ) 

EXPCON=DEXP(-6.8830220D+00*T/318 .80D+00) 

A(5) = _3 . 2559937 6D+ 00 +7 . 271143381D-01*T -2 . 04853587 6D+04*EXPCON 

A(4) =-7.3388518970+ 03 

A(3)= 1 . 284952625D+05 - 1 . 04054913D+02 *T +8 .78398364D+06 *EXPCON 
A( 2 ) =-l .064506759D+06 +1 . 170089157D+03*T -5.067990044D+07*EXPCON 
A( 1 ) = 8314 . 34D+00 *T 


DEXPCO= (-6.8830220D+00/318.80D+00) *EXPCON 

DA ( 5 ) = 7 . 271143381D-01 - 2 . 048535876D+04*DEXPCO 

DA (4) = 0.0D+0 

DA ( 3 ) =-l .04054913D+02 + 8 .78398364D+06 *DEXPCO 

DA(2) = 1.1700891 57 D+ 03 - 5 . 067990044D+07*DEXPCO 
DA( 1 ) = 8314 . 34D+00 


DDEXPC= 
DDA ( 5 ) = 
DDA { 4 ) = 
DDA ( 3 ) = 
DDA ( 2 ) = 
DDA { 1 ) = 
RETURN 
END 


( (-6 .8830220D+00/318 .80D+00) * *2. 0D+00 ) *EXPCON 

- 2 . 048535876D+04*DDEXPC 
0 .0D+0 

8.78398364D+06 *DDEXPC 

- 5 . 067990044D+07 *DDEXPC 
0. 0D+00 
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